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Abstract 

We study the O(N) vector model and the U(N) Gross-Neveu model with fixed total fermion 
number, in three dimensions. Using non-trivial polylogarithmic identities, we calculate the large- N 
renormalized free-energy density of these models, at their conformal points in a "slab" geometry 
with one finite dimension of length L. We comment on the possible implications of our results. 
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1 Introduction 



Conformed field theories (CFTs) in dimensions d > 2 have recently attracted much interest Q. Generic 
results based solely on conformal invariance are not very restrictive in d > 2, (see however |], |[]), 
and most of the information is extracted by studying explicit models. Although many CFTs have 
been recently discovered in d = 4 ||, a large amount of work has also been devoted to the study of 
CFT models in 2 < d < 4, such as the O(N) vector model || [?], |||] and the Gross-Neveu model ^ at 
their conformal points for large- N. 

In this work, we study the O(N) vector model and the U(N) Gross-Neveu model with fixed total 
fermion number, in d = 3. Using non-trivial poly logarithmic identities, we calculate the free-energy 
density of these models to leading order in the 1/N expansion, at their conformal points in a "slab" 
geometry with one finite dimension of length L. The free-energy density of the three-dimensional 



O(N) vector model in a "slab" geometry was first calculated in 11]. Recently, the free-energy 
density of (super)conformal field theories in d = 4 was also calculated in terms of poly logarithms [12|. 

Based on earlier ideas |13| , it has been recently argued in |l4]] that the free-energy density of four 
dimensional quantum field theories (QFTs) encodes non-perturbative information for the massless 
degrees of freedom coupled to a fixed point. Therefore, it may be interpreted as a measure of their 
expected reduction from the UV to the IR [p^[ . Our leading- N analytic results are consistent with an 
extension of such an interpretation for the free-energy density to three-dimensional QFTs. 

2 The Conformally Invariant O(N) Vector Model in d = 3 



We begin by reviewing the results of [11|. Consider the partition function of the O(N) vector model 



in d = 3, obtained after integrating out the fundamental scalar fields 4> a (x), a = 1, 2, .., N, 

Z B = J(Da)exp[-NS eff (a,g)} , (1) 

S eff (a,g) = lTr[ln(-d 2 + a)}-±Jd 3 xa(x), (2) 

where a(x) is an auxiliary scalar field and g is the coupling. Setting a(x) = M 2 + (i/y/N) cri(x), (||) 



can be calculated in a renormalisable 1/N expansion [16], provided the gap equation 
d 3 p 1 



g J (2tt) 3 P 2 + M 2 



(3) 



is satisfied. To any fixed order in 1/iV, a non-trivial CFT is obtained by tuning the coupling to the 
critical value 1/g = 1/g* = (27r)~ 3 / d 3 p/p 2 |lq]. Then, the renormalised mass (or inverse correlation 



2 



length) M(= l/£) = 0. 

When the model is placed in a "slab" geometry with one finite dimension of length L and periodic 
boundary conditions, the gap equation reads 

l l 00 r d? l 

with the momentum along the finite dimension taking the values uj n = 27rn/L, n = 0, ±1, ±2, ... We 
may then study the existence of a finite-temperature phase transition in the dimensionally continued 
version of the model, (i.e. in d — 1 infinite dimensions) fli~6| , |i~7| , |l8|] . Since UV renormalisation is 
insensitive to putting the system in a finite geometry, the coupling constant on the l.h.s. of (|j) can 
be set to its renormalised value in the bulk which explicitly depends on a mass scale Mo, such that 
the system is in its 0(A r )-ordered phase for zero "temperature" T ~ 1/L. Then, we can obtain an 
equation which gives the dependence of Ml on Mq and T. The finite-temperature phase transition 
corresponding to O(N) symmetry restoration, occurs when Ml = for some critical temperature T* 
which is related to Mq through the above equation. It can be shown that the finite-temperature phase 
transition cannot take place for 2 < d < 3 in accordance to the Mermin-Wagner-Coleman theorem, 



but can only occur for 3 < d < 4. The critical temperature T* obtained this way agrees [19] with 



the well-known results, i.e. see [2C]. The theory at T* is a d — 1-dimensional CFT, however the OPE 
structure of its correlation functions is less clear . 

On the other hand, when the coupling is fixed to its bulk critical value 1/g*, (|j) has the solution 

M L = M, = |ln(±±^ , (5) 

which corresponds to the physical situation of finite-size scaling |21] of the correlation length. Then, 
the subtracted Q free-energy density for this configuration reads 

foo-h 1 f d 3 p , 2 1 ~ f d 2 P , , 2 , 2 , , 1 f d 3 p Ml 



ln(e- LM *)Li 2 (e- LM *)-Li 3 (e- LM *; 



N 2 J (2tt) 3 "^ 2£ n f^ oo y (2tt) 2 "^ ' " n ' "~*' ' 2 J (2tt) 3 p 2 

M 3 1 



12tt 2ttL 3 

4 C(3) 

5 2?rL 3 ' 



(6) 



1 Since the UV singularities are the same in the bulk and in the finite geometry, subtraction of the bulk free-energy 
density ensures a UV- finite result. 
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where Li n (x) are the usual polylogarithms [22|. The third line in (|6|) follows from the second, by 



virtue of non-trivial polylogarithmic identities [p^j , and fits into the general formula [23] 



Joo - lL = C -g-jjg 5 (") 

for the finite-size scaling of the free-energy density in conformal field theories, with c/N = 4/5. It is 
quite remarkable that the value of c/N obtained from (||), turns out to be a rational number. This is 
reminiscent of central charge calculations in two-dimensional CFTs, (see |2f| for a recent reference). 
There is strong evidence f p6fl that correlation functions at the above finite-size critical point can by 
described by operator product expansions of the bulk O(N) vector model, in accordance with earlier 
ideas [23|. 



The result (^) is consistent with the interpretation of c as a measure of the massless degrees of 
freedom coupled to a critical point and their expected reduction from the UV to the IR. To see this, 
we recall that the large- iV critical point of the O(N) vector model is identical with the large- N IR 
critical point of the O(N) invariant 4> 4 theory p8|| . The latter theory has as UV critical limit the free 
theory of N massless bosons, for which it is well-known that c/N = 1 in d = 3. Therefore, according 



to the interpretation given in [14], we expect that 



cuv = c(N massless free bosons) > c(0(iV)vector model) = cir , (8) 

which is satisfied since 1 > 4/5. Note that (|8|) is a non-perturbative result from the point of view 
of the renormalisation group (RG) flow, since the large- iV IR critical point of the 4> 4 theory is non- 
perturbative in the coupling. 

3 The U(N) Gross-Neveu Model in d = 3 

Following the considerations above, it is possible to study the free-energy density of other CFT models. 
Consider, for example, the U(N) invariant Gross-Neveu model in d = 3 whose partition function, after 
integrating out the fundamental Dirac fermionic fields tp a (x), tp a (x), a = 1,2,.., N, reads 

Zf = jiVX) exp[-iV/ e// (A,G)] , (9) 
Ieff(X,G) = ^|d 3 xA 2 (x)-Tr[ln(^ + A)], (10) 



2 Recall that, in two dimensions c coincides with the central charge |24 
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where A(x) is an auxiliary scalar field and G is the coupling. We use the notation <f) = 7^9^ and 
the following two-dimensional Hermitian representation for the Euclidean gamma matrices in d = 3 

I27LEI 



71 = a ! 72 = cr , 73 = 70 = cr , (11) 

where a*, i = 1,2,3 are the usual Pauli matrices. This model describes fermion mass generation 
through the breaking of space parity |e3] . The partition function (El) can be evaluated in a renormalis- 



able 1/N expansion 27] when one sets X(x) = m + (1/s/N) Ai(x), provided the following gap equation 
is satisfied 

1 - 2 f d3p 1 (12) 
G J (2tt) 3 p 2 + m 2 ' 1 ; 

At the critical coupling 1/G = 1/G* = 2(2ir)~ d f d d p/p 2 , the theory is conformally invariant and 
m = 0. 

When the system is placed in a "slab" geometry with one finite dimension of length L, the fermions 
acquire antiperiodic boundary conditions and the gap equation reads 

i = 2 " r d^p i 

G L n yj (27r)2p2 + a ;2 +m | > K > 

with oj n = (2n + 1)tt/L, n = 0, ±1, ±2, ... Again, we may study the finite-temperature phase transition 
in the dimensionally continued version of the model, (i.e. in d — 1 infinite dimensions) |]27| , O ], by 
setting 1/G to its bulk renormalised value which explicitly depends on the mass mo of the fundamental 
fermionic fields. This means that the system is in its "broken phase" for zero "temperature" T ~ 1/L. 
Then, we can obtain an equation which gives the dependence of mi on mo and T. The second order 
finite-temperature phase transition corresponding to space parity restoration, occurs when = 
for some critical temperature T* which is related to mo through the above equation. The finite- 
temperature phase transition is now possible for all 2 < d < 4 p7| , O], due to the absence of zero 
modes for fermions and antiperiodic boundary conditions. 



On the other hand, when the coupling stays at its bulk critical value 1/G*, (13) is satisfied for 
rriL = m* = 0. (14) 
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This essentially means that, to leading order in 1/N, the free-energy density of the system is given by 
the free-field theory result.^ Indeed we easily find 

foo - h _ 3 C(3) n -. 



which implies that c/N = 3/2 in agreement with the results of [29]. 

Although the result (|l5|) may seem trivial, it encodes important information. This is again seen if 
we recall that the large- AT critical point of the U(N) Gross-Neveu model is identical with the large- N 
IR critical point of the Gross-Neveu- Yukawa model p8[ . The latter model has as UV critical point 
the free theory of N massless Dirac fermions plus one massless boson for which c/N = 3/2 + 1/N in 



d = 3. Then, according to [14[ we expect 



cuv = c(N massless Dirac fermions + 1 massless boson) > c(Gross-Neveu) = cm , (16) 



which is consistent with (|i5|) for large and finite N. Again, the non-perturbative nature of ( |i6[ ) is a 
consequence of the non-perturbative nature of the large- IR critical point of the Gross-Neveu- Yukawa 
model. 

4 The U(N) Gross-Neveu Model in d = 3 with fixed total fermion 
number 



The Gross-Neveu model can also be studied p0[ |31[| for fixed total fermion number B. To this effect, 

we introduce a delta-function constraint 5(N — B) into the functional integral @ at finite temperature, 
where 

N = J d 2 xV f (x)V(x), x=(xi,x 2 ) (17) 

is the fermion number operator. Using an auxiliary scalar field 9{xs), the above delta-function con- 
straint is exponentiated and after integrating out the fermions we obtain the partition function 

Z f = J(pX)(pe)esxp[-NX e f f (X,g;0,B)] , (18) 
T e// (A,g ; 0,S) = i S^(x3)d 2; 3 + ^^d 3 xA 2 (x)-Tr[ln^ + i73^ + A)] L , (19) 

where B = B/N is assumed to be finite for large- N and \/Q is the new coupling. The subscript L 
denotes X3-integration up to L and the latter quantity plays here the role of inverse temperature 1/T. 
3 There exist solutions of ( |l3| ) with imaginary toil which will be discussed elsewhere. 
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For large-TV, the functional integral (18) can be calculated by the steepest descent method, since 
it is dominated by the uniform stationary points (A) and (6) of I e ff- These stationary points are 
obtained as the solutions of the following set of saddle-point equations 



1 

G 



r> OO 

tE 



d 2 P 



LJ^J (2vr) 2 p 2 + {u n + (6)) 2 + (A) 2 ' 



(20) 



ib 



lim — 

r^O L 



E 



d 2 P 



(2tt) 2 p 2 + K + <0» 2 + (A) 2 



(21) 



where b = (BL/tt), with f2 the total volume. The regulating term e lu3nT on the r.h.s of (|2l| ) has been 
discussed in |3lJ and ensures a finite result in the limit r —* 0, after the Matsubara sum has been 
performed. For (9) purely imaginary, which corresponds to having a real chemical potential \x = —i{6) 
p0| , j3lf| , (P0|) coincides with a similar saddle-point equation obtained in |l7], Wa . We can renormalise 
(p0|) by substituting for 1/Q the bulk renormalised coupling 1/G from (|l2|), since the presence of (9) 
does not alter its UV behavior. In this way, we can study the finite-temperature phase transition of 
the model in terms of the renormalised mass of the bulk fermionic fields and the chemical potential. 



For example, (21) would now give the critical fermion number B cr at which space parity is restored. 

However, we are interested here in possible real values of (9) which satisfy ([H]) and (^). The reason 
is that, if (9) is a real number we can set in (p0|) the coupling 1/Q to its bulk critical value 1/G* and 
obtain the following equation for (A) 

L (A) + In (l + e-W-WQ) + In f 1 + e~ L ^ +iL ^) = . (22) 



This has a real solution for (A) in terms of L, whenever we have —1 < cos (L(0)) < —1/2, or simply 
2vr/3 < L(9) < ir. Note now that (A), which is the renormalised inverse correlation length, is non-zero 
for L(9) 7^ 27r/3 corresponding to a finite-size scaling regime for our fermionic model. Such a regime 
is absent for the usual Gross-Neveu model studied in Section 3 and it is a consequence of keeping 
the total fermion number fixed. The dimensionless quantity L{\) is plotted in Fig. 1 for the allowed 
values of (9). 



The second saddle-point equation (21) ensures a fixed mean fermion number (N) = B, as imposed 
by the constraint. It turns out that its r.h.s. is real and has a quadratic divergence which should be 
subtracted. Again, if (9) is purely imaginary ( ^T]) assumes the usual form |32|] of the expression for the 
conserved charge in a system exchanging particles with a reservoir. It also vanishes for fi = —i{6) = 0, 
as it should, since this would correspond to the absence of conserved charges. On the other hand, for 
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7 n 



<e> 



Figure 1: L(X) for the allowed region of (6). 



2ir/3 < L(9) < ir we obtain, after subtraction of the divergence, 



2vrL 2 



[Cl 2 (2<j)) - Cl 2 (2cj) - 2L{9)) - Cl 2 {2L{9))] 



arctan 



e -L(X) gin (£-(#)) 
1 + e -L(A) cos (L(9)) 



(23) 
(24) 



where ClzioS) = Im [Li2(e ltJ )] is Clausen's function p5]. Now, b is related to the total fermion number 
of the system and in principle it should be real and positive. Therefore, from (p3|), the only allowed 



real value for (9) which satisfies both saddle-point equations (^y) and (21) is (6) = tt/L. In this case, 

6 = and it seems that there are no fermions left in the system. This is consistent with the apparent 
bosonization of the theory for (9) = n/L which will be discussed shortly. However, it is conceivable 

that the imaginary solutions for b may also have physical meaning, as they give rise to a real value for 
the free-energy density of the theory. The latter result is rather surprising and is obtained by virtue 
of non-trivial poly logarithmic identities. 

To demonstrate the above points, we calculate the free-energy density of the model and we obtain 



foo jh 
N 



d 2 p 

/. — ./ (2tt) 2 



-i oo 

7 E 



ln(p 2 + K + <0» 2 + (A) 2 ) 



d 3 p 2 (A) 2 



(2ny 



■ lnp 



2G* 



i(6)b 



In 



+ - 



d 3 p 
(2tt]3 

1 f d 2 p 



p 2 + (A) 2 \ (A) 2 



LJ (2vr) 2 
d 2 p 



-i{9) 



l n j 1 + e -LVp 2 +<A> 2 -iZ<f)j + l n ^1 + e -Wp2+{\)2+iL(9) 
1 1 1 



(2vr) 2 \ 1 + e W P 2 +(\) 2 -iL(e) 



1 + e 



L v / p 2 + (A> 2 + l L(e) 
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(A) 3 + 1 



+ 



6-7T 7T L 3 

HO) 



lnfe^W) Li 2 (-e-^>,L(0)) - Li, (-c~ L W,L(9)) 



2vrL 3 



C7 2 (2<£) - CZ 2 (20 - 2L(0)) - Cl 2 {2L(9))} . 



(25) 



Li n (r, 9) is the real part of the polylogarithm Li n \ re l ®j in Lewin's notation [p^ j. As already mentioned, 
the free-energy density ( ^5| ) is real for 27r/3 < < vr. 

For (0) = tt/L, (|2^ ) has the solution (|5j) (A) = (the "golden mean") which corresponds to the 
physical situation of finite-size scaling of the correlation length in the O(N) vector model in d = 3. 
This apparent "bosonization" of the Gross-Neveu model is effectively a transmutation between Fermi 
and Bose statistics, since the solution (6) = tt/L introduces a zero mode for the fermions. Also, for 
this particular value of {9), b is zero which can be further interpreted as the absence of conserved 
charges in the bosonized version of the system. The bosonization of the model for (9) = tt/L is also 



discussed in 33 1 and a possible connection with anyonic physics is made in [34]. In this case, using 
the same polylogarithmic identities as in @ we obtain 



foo Jl 

N 



8 C(3) 
5 2vr L 3 ' 



(26) 



which is consistent with the expected CFT result ([?]), with c/N = —8/5. In fact, the full expression 
(pq) is consistent with the scaling form ([?]) for the allowed values of (9), however the corresponding 
expressions for c/N involve polylogarithms and Clausen's functions at non-exceptional arguments 
and are not illuminating. In Fig. 2 we plot the numerical values of c/N corresponding to ( |25|) for 
2vr/3 < L(9) < vr. 



7 J[ 



<e> 



Figure 2: c/N for the allowed region of (9). 

Nevertheless, at the other end-point of the allowed (9) region, (9) = 2ir/3L, (p^) has the solution 
(A) = which corresponds to infinite correlation length. Therefore, we may associate this point with 
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a temperature phase transition. Remarkably, (25) simplifies considerably for (0) = 2tt/3L giving 



foo fh 



N 2vrL 3 



(27) 



It may be interesting to point out that Cliiji /3) is the absolute maximum of Clausen's function [p2| 
which is a well-documented numerical constant. 

Consistency of our main result (25) with the expected from conformal invariance scaling form (0) 



raises the possibility that, at the above finite-size critical point with 2n/3 < L{0) < tt the Gross-Neveu 
model with fixed total fermion number is related to conformal field theory. This is most clearly seen 
for (0) = tt/L, when from ( p6| ) one concludes that the model is related to the O(N) vector model at 
its finite-size scaling critical point. Note that the free-energy density given by (^6|) is negative which 
shows that such a critical point is unstable by itself. 

The interpretation of our results (^) and (|27j) in terms of the loss of degrees of freedom from the 
UV to IR is not clear, except perhaps at the point {0) = tt/L which is related to the critical point of 
the O(N) vector model in d = 3. In particular, the numerical value of c in ( |27| ) is ~ 2.9 and it is larger 
than the sum of c's at the UV fixed points of the (j) theory and the Gross-Neveu- Yukawa theory, when 
N > 2.5. The latter sum may be considered as an upper bound for c in view of the known universality 
classes in three-dimensions. Therefore, the point {0} = 2ir/3L corresponds to a genuinely new phase 
transition for which the inequalities suggested in |14| do not seem to hold. Such a critical point may 
be associated with a non-unitary CFT and further with a Lee- Yang edge singularity p5|]. In fact, an 
interpretation of i(0) as a purely imaginary chemical potential suggests |3(| that, for 2-7r/3 < L(0) < tt, 



we are dealing with critical points corresponding to Lee- Yang zeros [37]. 



5 Summary and Outlook 

In this work we calculated the renormalised free-energy densities of the three-dimensional O(N) vector 
model and the U(N) Gross-Neveu model, at their conformal points in a geometry with one finite 
dimension of length L. Our calculations were based on non-trivial polylogarithmic identities. Our 



results are consistent with the recent suggestion of \14] that the free-energy density encodes non- 
perturbative information regarding the loss of degrees of freedom from the UV to the IR in QFTs. 

Our study demonstrates the importance of three-dimensional critical models as testing grounds for 
the recent ideas concerning irreversibility of the RG flow and non-perturbative phenomena in QFTs 
|j~4| , pq| . There are many directions in which our leading- iV calculations could be extended. For 
example, one could study the IR limits of the O(N) vector model and the Gross-Neveu model in 
three-dimensions, both in the bulk and in a finite geometry. These limits would correspond to the 
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lower bounds of the inequalities suggested for c. The latter could be further checked by next-to-leading 
order calculations of c in the above models^]. Based on techniques developed here, one could also study 
the free-energy density of other three-dimensional models exhibiting critical behavior, such as CP N ^ 
models @ or QED 3 Q. 

In an possible application of our results, we may view the critical points for 2ir/3 < L(9) < ir as viable 
critical points in supersymmetric CFTs in d = 3. For example, consider the N = 1 supersymmetric 
cj-model in d = 3 [42] which contains Majorana fermions. This means that the large-TV fermionic 
contribution to the free-energy density of such a model with fixed total fermion number is half the 
r.h.s. of (^), therefore the sum of the bosonic and fermionic contributions ^3|, as seen from (0) and 
Fig. 2, is always greater than zero. This sum vanishes for (9) = tt/L corresponding to a supersymmetry 
restoration mechanism which deserves further study. 

The relevance of non-trivial polylogarithmic identities to the calculation of free-energy densities in 



two-dimensional CFTs is well-known [25]. Their appearance in studies of CFTs in higher dimensions is 
intriguing and requires further investigation. It would also be interesting to study the OPE structure 
of correlation functions in the Gross-Neveu model for fixed total fermion number, at the above critical 
points. Although our results simplify for d = 3, they can presumably be generalised for all 2 < d < 4. 
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